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We study adiabatic charge pumping through a quantum dot placed at the junction of N quantum 
wires. We explicitly map out the pattern of pumped charge as a function of the time-varying 
tunneling parameters coupling the wires to the dot and the phase between any two time varying 
parameters controlling the shape of the dot. We find that with N — 2 time-independent well-coupled 
leads, the maximum pumped charge in the remaining two leads is strongly suppressed with increasing 
N, leading to the possibility of tuning of the pumped charge, by modulating the coupling of the 
TV - 2 leads. 

PACS numbers: 73.23Hk,72.10Bg,73.40Ei 
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An adiabatic charge pump[l| was theoretically pro- 
posed by Thouless for an isolated system and later ex- 
tended to open systems at finite temperature @, 0, lij. It 
is essentially a device that drives current through a sys- 
tem at zero bias by a time-dependent variation of two or 
more parameters. Although theoretically proposed about 
a decade ago, Switkes et al\j^ recently demonstrated that 
an adiabatic electron pump, which produces a dc voltage 
in response to a cyclic change in the confining potential of 
an open quantum dot was experimentally feasible. This 
has led to renewed interest @, 0, H, @, E U, HI, IU 14, 
H, in the field. 

In this paper, we study adiabatic charge pumping 
through a quantum dot placed at a junction of N quan- 
tum wires and which is well-coupled to the wires. We 
show that under a time dependent variation of two of the 
tunneling parameters to the dot, charge Qij = Q(i,r) — 
Q(j, t) can be pumped from wire i to wire j. An impor- 
tant point to note is that unlike the case when N = 2, 
the transmission between any two wires for N > 3 never 
reaches unity anywhere in the parameter space, because 
of the possibility of leaking out through the other wires. 
Thus, analogously to the N = 2 caseflij. large pump- 
ing occurs when the 'pumping contour' (a closed curve 
in the Ji-Jj plane, where Jj and Jj denote the tunneling 
parameters between wire i,j and the dot) encloses trans- 
mission maxima peaks (when Y] ■ — 1) and cuts the 
resonance line at points where the transmission is small, 
thus ensuring that most of the resonance line is enclosed. 

Consider a time-dependent scatterer connected to N 
single channel leads. In the adiabatic limit, at lead- 
ing order, (O(w ), where u> is the frequency), the in- 
coming particle sees a static scatterer, and hence the 
scattering and the current can be computed using time- 
independent scattering theory (the Landauer-Buttikcr 
formula). But as shown by Buttiker, Pretre and 
Thomas [J] and BrouwerQ, even at order O^ 1 ), the cur- 
rent and the charge pumped can be computed from the 
'frozen' or 'quasi-static' scattering matrix[l2] data. For 
a simultaneous slow and small variation of two time- 



dependent parameters X\ and X2, 
from wire i is given by 
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where the first term is denoted as Qmt and the second 
as Qtran and the S'-matrix for an A^-lead system is given 
by 
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Ta , tij , 8i and ipij are time dependent in general. We will 
mainly use Eq.@ in our computations. Note that the 
first term on the RHS of Eq. ^ is a total derivative and 
denotes the winding number around the R = singular- 
ity. Thus Qi n t is forced to be an integer. The second 
term is not quantised, however, and hence, the pumped 
charge is not an integer. But by choosing the contour 
appropriately, to cut the transmission resonance maxima 
when it is small, Qtran can be made small, and the charge 
pumped is almost quantised [l4|. As we shall see, this is 
harder to arrange when N > 2, because there is a possi- 
bility of transmission into the other wires. This explains 
why the value of the maximum pumped charge decreases 
as N increases, which is one of the main results of this 
report. 

Let us now write down the lattice model for a quantum 
dot at the junction of N quantum wires. The wires are 
modeled as a one-dimensional chain of sites, with on-site 
energy zero and with nearest neighbour hopping ampli- 
tudes given by — J. The dot is again modeled as a bunch 
of tight- binding sites with nearest neighbour hoppings 
given by —Jo and on-site energy eo. Finally, the dot is 
coupled to the i th lead by a hopping element J^. In the 
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work presented in this brief report, we have taken Jjj = J 
and e = for simplicity, since taking Jn ^ J and £ ^ 
makes no qualitative difference to the resultsfTij]. Thus, 
the Hamiltonian of the system is given by 

H -Pwircs ~t~ E^dot — wire ~t~ -^dot 

N L 

= [ ~ J C lrn C irn+l ~ ^4^+1 

i— 1 m—l 
l-l 

-J Y c tm c *m+i - J4>c a + h.c] (4) 

m—l 

Here L denotes the number of sites on each wire with 
wire index i and I denotes the number of sites on each 
wire representing the dot. cq denotes the dot point which 
is common to all the wires. Since our focus here is on 
charge pumping, we have not included spin in the above 
Hamiltonian. 

For fixed values of J, and Ji, (i = 1, . . .N), the scatter- 
ing problem is straightforward. Let us assume that the 
incoming wave is along wire 1 , without loss of generality. 
Then the wave-function at the site n in lead 1 is given by 
Ci n = cxp(—ikx n ) + 7*11 exp(ikx n ); ( n > I + 1) and the 
wave- function at site n in all other leads i = 2, . . . , N 
are given by Ci n — tu exp(ikx n ); ( n > I + 1) where 
x n = na for the lattice spacing a. Note that the num- 
bering of the sites, implies that incoming waves are de- 
noted by exp(— ikx n ) whereas outgoing states are de- 
noted by exp(ikx n ). The wave-function at sites within 
the dot is given by (j = 1, . . . N) a jn = P( e^q>(—iqx n ) + 
P% exp(iqx n ); ( n < I). The wave- vectors k and q are re- 
lated to the energy as 

Ek = — 2Jcosfca and Ek — £o — %Jd cosqa (5) 

as can be seen by solving the equations in the bulk of 
the wire and the dot. The set of boundary equations (for 
amplitudes at the boundary of the dot and the wire, i.e., 
sites / and I + 1 and at the junction point of the dot) can 
be solved to obtain r\\ and tu in terms of J, Jd and Ji, 
but it gets progressively more difficult for higher values of 
N to get analytical results. However, the equations can 
be solved numerically and we can obtain the scattering S- 
matrix in the parameter space of the hopping amplitudes. 

To obtain pumped charge, any two of the hopping pa- 
rameters that couple the dot to the wire need to be adi- 
abatically varied. Hence, in the lattice model, two of the 
hopping parameters that couple the dot to the wires are 
taken to oscillate with the frequency to, with a modu- 
lation parameter Pi and a phase difference 2<f>, i.e., we 
choose 

Ji = Jm + Pi cos(ujt + <t>) 
J =ti = Jjo + Pj cos(ujt - 4>) , (6) 

where j is one value not equal to i. In this report, we 
choose Pi = Pj. The case where Pi ^ Pj, as also the 



cases when more than two parameters are varied will be 
considered elsewhere fl8| . 

The charge pumped out from wire i, Q(i,r) is then 
computed from the time dependence of the S'-matrix pa- 
rameters ru, which, in turn, are known in terms of the 
hopping parameters, using EqO Finally, we can obtain 
the charge pumped between different wires as 

Q ij = Q{i,T)-Q(j,T) . (7) 

Note that Q(i,r) is computed with the wave incident 
in wire i, whereas Q{j,T) is computed with the wave 
incident in wire j. The results for the pumped charges 
for the different cases are given below. 

The parameters that affect the value of the pumped 
charge are the pumping amplitude P and the phase <fi. 
We have varied <j> over its full range. Within the adiabatic 
approximation, the pumped charge is independent of u) 
and so we have chosen a specific value of ui. 

We have solved the set of boundary equations for N = 
2,3 and 4 wires. For N = 2, we have chosen IN = 4, 
and both the hopping amplitudes coupling the dot to the 
wire are made time-dependent using Eq(6j and we obtain 
a single pumped charge Q = AQ12 — AQ2i- For N = 3, 
any two of the three hopping amplitudes coupling the dot 
(with IN +1 = 4 sites) to the wire can be made time- 
dependent. Let us call them J\ and J2. We then obtain 
the pumped charge Q — AQ12 — AQ2i- (Similarly, by 
making the other hopping amplitudes time dependent, 
we can obtain Q 1 — AQ23 — AQ32 and Q" = AQ31 — 
AQ13). The pumped charge Q is computed when the 
third hopping amplitude J3 = l.f which is of the same 
order as the mean values of Ji^ within each cycle. 

We have plotted the pumped charge Q as a function 
of the pumping strength Pi. 2 = P and the phase (f> m 
the contour plots of FigQ] for TV = 2 and 3. We have 
chosen (j> to go from — * ir, since EqJS] shows that the 
difference in the phase between the tunneling parameters 
on the two leads connected to the dot, is given by 2</>, 
and we have chosen a maximum of 10 for P. In Fig (21 we 
have plotted the different classes of contours in the P — <p 
plane for N = 2 and 3. We find that the contours can 
be classified into 3 catagories. Those which (A) do not 
enclose either of the two |T| = 1 points, (B) includes one 
of them more number of times than the other and finally 
(not shown for TV = 3) (C) includes both of them equal 
number of times. We see (as was shown in Ref. [141 ]) 
that the pumped charge is large only in the case (B) 
with the sign of the pumped charge depending on the 
sense in which the contour encloses the \T\ = 1 point. 
We note that the pattern of pumped charge remains the 
same for N = 2 and 3, since the form of the variation 
of the paramters remains the same. But the magnitude 
of the pumped charge decreases for N = 3. This is due 
to the 'leak' of the wavefunction in the third direction, 
so that the transmission never reaches unity in the 1-2 
direction. 
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FIG. 1: Contour plot of the pumped charge for N = 2 and 3 wires going from left to right. The x- and y- axes are Pcos<f> and 
Psin <j>. The colour scale is on the right of each plot. Note that the pattern remains the same in both ; for TV = 2, the pumped 
charge reaches a maximum of ~ ±2, whereas for N — 3, it reaches a maximum of ±1.11. 
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FIG. 2: Typical pumping contours as a function of the pump- 
ing parameters Xi = J? and X2 = jf- For N — 2, the full 
pumping contours are shown for three classes corresponding 
to the regions A,B and C in Fig. 1(a). For N = 3, only a 
small part of contours of classes A and C corresponding to 
regions A and C in Fig. 1(b) is shown, along with the reso- 
nance points R = 0. Note that contour A does not enclose 
the resonance points and contour C encloses both once 

In Fig. 3, we have separately shown the contributions 
of the the winding number term Qi n t and the transmis- 
sion term Qtran(see Eq.Q for N — 2 and 3. These 
graphs clearly show that large pumped charge is always 
correlated with the winding number term being non-zero. 
Fig. 3(b) also shows that the range of values of P for 
which the pumped charge is large is narrower for N = 3. 
This can be correlated with the coming together of the 
resonance lines as the coupling to the third wire is in- 
creased. This is shown explicitly in Fig. 4(a), and also in 
Fig. 2(b), where it is clear that it is hard to find contours 
which enclose one resonance but not the other. Thus, the 
parameter range (P c , 4>), for which the contours enclose 
one maxima and not the other is quite small, and so the 
region where pumped charge is large is small. 

In FigfJJb), we have shown the plots of the pumped 
charge when J3 <C J\,i and J3 3> Ji.2- In either case, 



FIG. 3: The different contributions to the pumped charge Q 
for N = 2 (a) and N = 3 (b). The full red line is Q int , the 
blue dash-dot line is the total pumped charge Q and the green 
(and magenta for TV = 3) dotted lines are Qtran- Note that 
large pumped charge is correlated with Qint 7^ 0. In Fig. 3(b), 
only a small region of P is shown, because the width of P for 
which pumped charge is large is smaller. The numerical values 
chosen to obtain Q are Jd — J — 1.0, J10 = J20 = 2.0,0 = 
0.057rfca = O.OOItt, u> = 1.0. 

it can be checked that the plots are very similar to the 
N = 2 case studied earlier This is not surprising 

because both for J3 very small and J3 very large, the 
leak to the third wire is very small. But for intermediate 
J3, we see that not only is the maximum pumped charge 
smaller, it occurs for a narrower range of parameters as 
expected. 

What happens when 4 wires are attached to the dot 
and two of the hopping amplitudes coupling the wire 
to the dot are made time-dependent? Can we identify 
trends as the number of wires N increases? To check this, 
we obtained the numerical results for N = 4 wires and 
we find that the magnitude of pumped charge decreases, 
with regions of large pumped charge getting thinner as 
we increase N. In FigfSl we plot the maximum pumped 
charge as a function of the number of wires well-coupled 
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FIG. 4: (a) Contour plot of transmission maxima (T12 + T13) 
for three values of J3. Lines A, B and C denote resonance 
lines for J3 = 0.1.0.9 and 1.1 respectively. The crosses denote 
the R = points. Note that in the lower left quandrant, lines 
A and B are indistinguishable and the resonance point on 
A,B and C is indistinguishable and denoted by a single cross, 
(b) Pumped charge Q as a function of the pumping amplitude 
P for three values of J3 = 0.1, 1.1 and J3 = 10.0. The other 
parameters are the same as in Fig[2] 

(and held constant in time) to the dot. 

Hence, we conclude that it is harder to tune for 'almost 
quantised' pumped charge for pumping through a dot 
at a junction of N wires except for N = 2. However, 
the magnitude of pumped charge is not vanishingly small 
for N = 3,4. There is a suppression of the maximum 
pumped charge from ~ 2 for N = 2 to about 1.2 for N — 
3 and to 0.45 for N = 4. This can be easily explained 
as a 'leak' of the electrons into the wires that are not 
participating in the pumping process. Hence, it is easy 
to see that no pumped charge is expected in the large N 
limit. 

In conclusion, we have explicitly obtained the pumped 
charge as a function of the various parameters such as 
pumping amplitude, phase <fi, pumping frequency, etc, for 
pumping through a quantum dot placed at the junction 
of N wires where TV = 2, 3 and 4. We note that the maxi- 
mum pumped charge reduces in magnitude as the number 
of wires is increased because some of the charge can leak 
out through the other wires which are not involved in the 
pumping process. The pumped charge is also seen to be 
large only when the pumping contour encloses the trans- 



mission maxima (not necessarily a \T\ = 1 point) appro- 
priately and cuts the resonance line at a point where the 
transmission is low. 

More general cases which are possible for N > 3, 
which involve time dependent variation of all the tun- 
neling parameters coupling the wire to the dot can also 
be studied [l^| . One might expect that unless the time- 
dependences are carefully adjusted the pumped charge 
would be quite low. However, they have not been explic- 
itly considered here. The extensions to higher uj, where 
the adiabatic approximation breaks down will also be 
presented elsewhere. 
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